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Abstract. The expected episodic or intermittent star formation histories (SFHs) of gas rich dwarf irregular galaxies (dlrrs) 
are the longstanding puzzles to understand their whole evolutional history. Solving this puzzle, we should grasp what physical 
mechanism causes the quiescent phase of star formation under the very gas rich condition after the first starburst phase. We 
consider that this quiescent phase is kept by lack of H2, which can be important coolant to generate the next generation of stars 
in the low-metal environment like dlrrs. Furthermore, in dlrrs, H 2 formation through gas-phase reactions may dominate the 
one on dust-grain surfaces because their interstellar medium (ISM) are very plentiful and the typical dust-to-gas ratio of dlrrs 
(Odirrs = 1.31 x 10 -2 Dmw, where £>mw is its value for the local ISM) is on the same order with a critical value D a ~ 10~ 2 £>mw- 
We show that the lack of H2 is mainly led by H destruction when gas-phase H2 formation dominates since H is important 
intermediary of gas-phase H2 formation. H is destroyed by the radiation from all stars born in the previous starburst phase 
because Hr destroying infrared photon can penetrate the whole ISM of dlrrs. Considering the physical process which timescale 
is the longest as main process in regulating global star formation, we can show this lack of H2 leads the quiescent phase of 
star formation. Hence, we can say that the stellar radiation which destroys H and leads low H2 abundance should be properly 
treated in studying SFHs of dlrrs. 
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1. INTRODUCTION 

Dwarf galaxies are the galaxies with lower luminosity, lower 
mass and smaller size, and they are much more abundant in 
number than normal galaxies on the Hubble sequence like our 
Milky Way (see Kunth & Ostlin 2000, for a recent review of 
local group dwarf galaxies). These natures imply that they are 
basic building blocks of the normal galaxies in hierarchical 
structure formation scenarios of the cold dark matter u niverse 
jCole et alJl2000HCorbin & Vaccal20o3lMonical2003l) . In ad- 
dition, the dwarf galaxies have lower abundances of heavy ele- 
ments than those of the normal galaxies; therefore, they are still 
in an early stage of che mical evolution and/or a re analogies to 
primordial galaxies (e.g.. llzotov & Thuanll999l) . 

Dwarf galaxies are classified by their morphological ap- 
pearances into some basic types. One of the morphologi- 
cal types is dwarf irregular galaxies (dlrrs), which charac- 
teristics are irregular and amorphous appearances at optical 
wavelengths. Interestingly, unlike other morphological types 
of dwarf galaxies, dlrrs still have plentiful interstellar medium 
(ISM ) and in general show ongoing star f ormation and HII re- 
gion (Mateo 1998; Kunth & Ostlrnl l2000l) . The star formation 
histories (SFHs) of dlrrs are considered to be episodic, that 
is, quiescent phase follows a starburst phase jGreggio et alJ 
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119931 lTolstovetalJll998t lOrebell Il999t) . and their star for- 
mation rates (SFRs) are more widely distributed (~ 10~ 4 - 
1 0~ 2 [Mfri / yr /kpc 2 l ) than that of the normal galaxies jHunterl 
119971) . In this issue, we particularly focus on the dlrrs which 
are considered to be in their quiescent phase of star formation 
(i.e. SFR ~ 10~ 4 [M /yr/kpc 2 ]). 

However, despite their closeness, the regulation of star for- 
mation even in the Local Group dlrrs is not yet well under- 
stood. Actually, it is very surprising that they have been form- 
ing s tars recently because their gravitational potential are shal- 
low (Skillman & Bender 1997, for a review) and comparable 
to the energies of some supernova explosions; therefore, they 
should have little ISM after some supernovae (several million 
years from the epoch of dlrrs formed), and thus, the recent and 
ongoing star formation would not occur. Why do they still have 
plentiful gas? Why do they show ongoing star formation? Both 
the questions are still controversial. Even if these questions are 
solved, another puzzle remains; why do they form stars at such 
low rates in spite of having plentiful ISM as material of stars? 
This paper attempts to solve this puzzle by investigating the 
non-chemical equilibrium molecule formation history of dlrrs 
on some assumptions. 

For the star formation processes in low-metallicity environ- 
ment, hydrogen molecules (H?) can be important cool ants at 
the gas temperature T $ 10 4 K JPeebles&Dickelll 968); more 
abundant H2 leads gas to cool faster, and therefore, would re- 



MARK & H. Kamaya: The self regulating star formation of gas rich dwarfs 



2 

suit in more active star formation. Although H2 forms primarily 
on the surfaces of dust grains in local molecular clouds, H2 can 
also form in gas-phase through the reactions, which intermedi- 
aries are negatively charged hydrogen (H ), 

H + e -> H + y, (1) 
H + H -> H 2 + e , (2) 

or ionized hydrogen molecule (HJ), 

H + H + -> H+ + y, (3) 
H++H -> H 2 + H + . (4) 

These gas-phase reactions are active T ~ 10 4 K and the 
reactions Q and can more efficiently form H2 than the 
Q and © reactions (e.g., lAbel et all Il997t) . On the other 
hand, H2 is dissociated by collisions with other chemical 
species and the absorption of ultraviolet (UV) photons called 
the Lyman-Werner band (LW) photons (11.26 eV < hv < 
13.6 eV), which are mainly radi ated by early type stars 
(Field. Somerville & Dressier 1966). 

Many astrophysicists have built the theoretical models con- 
centrated on the dissociation of H2 by the LW photons because 
they have investigated the second epoch of star formation after 
the fi rst stars have been formed in met al-free molecular cloud 
(e.g., Isiii3ll977t iNishi & TashirolbOOOll . These models result 
in the low abundance of H2 in relatively large region (several 
hundred parsec around the source star of the LW photons), and 
thus, negative-feedback on further star formation in the host 
cloud of the star. However, in dlrrs, the region affected by LW 
photons would be smaller than previously predicted because 
dlrrs are gas rich system; the mean free paths of LW photons 
are on the order of 1 pc (<§: the typ ical size of dlrrs) with typica l 
parameters of primordial clouds (Hai man. Rees & Lo eb 1996). 
Furthermore, early type stars as the sources of LW photons are 
rare in dlrrs. Therefore, the negative-feedback on further star 
formation brought about the H2 dissociation by the LW pho- 
tons may be negligible in dlrrs. 

By the way, if H2 is formed by gas-phase reactions (equa- 
tions ! 1I4> . the formative abundance of H2, yn 2 , largely depends 
on the H abundance, yu-, because the H2 formation process 
through H is more effective than the one through HJ (in §2, 
we discuss that gas-phase H2 formation may dominate grain- 
catalyzed H2 formation in dlrrs). We expect that yn 2 in dlrrs is 
more effectively influenced by the destruction of FL than by 
the direct dissociation of H2 due to the LW photons. This is be- 
cause that the typical mean free path of the H destroying in- 
frared (IR) photons (its wavelength A ^ 1 .64 fim), Zir, would be 
sufficiently longer than that of the LW photons; Zir ~ 100 kpc 
(» the typical size of dlrrs) with typical parameters of primor- 
dial clouds. Therefore, the IR photons can penetrate and affect 
the whole ISM. Furthermore, the IR photons are also radiated 
by low-mass stars as well as early type stars. Thus, all stars born 
in the previous starburst phase can be considered as the sources 
of the H destroying photons. If FF is seriously destroyed by 
the IR photons of the interstellar radiation field (ISRF) of dlrrs, 
the star formation activity may become more inactive at almost 
the whole ISM. 



With this expectation, we construct the theoretical model 
to investigate y^. In §2, we compare the rates of gas- 
phase and grain-catalyzed H2 formation, and show gas-phase 
H2 formation may dominate in dlrrs. §3 describes our model 
for the star formation process in dlrrs and §4 presents our re- 
sults. Finally, we summarize our study and briefly discuss about 
the next starburst phase in §5. 

2. GAS-PHASE AND GRAIN-CATALYZED 
H 2 FORMATION 

H2 formation on dust-grain surfaces completely dominates gas- 
phase one in the local ISM. However, if the dust-to-gas ratio, D, 
is lower than a critical value, D CI , g as-phase H? formation dom- 
inates grain-c atalyzed formation dKamava & Hirashitall200lt 
Kjloverl 120031) . Despite the H> forma tion rate on grain sur- 
faces is still uncertain, Glover (2003) calculates D CI in some 
astrophys ical situations by tentatively adopti ng t he rat e of 
iHollenbach & McKed Jl979l) . As noted in iGloverl d2003l) . the 
temperature dependence of D a is very strong; at low temper- 
atures (T £ a few x 10 2 K), grain-catalyzed H2 formation is 
relatively efficient and D a is very small (D cr S 10~ 3 Dmw, 
where £>mw is the dust-to-gas ratio in the local ISM), while 
at high temperatures, gas-phase H2 formation is efficient and 
thus, D cr is large. In addition, D cl also depends on the ionisa- 
tion degree and the density of ISM; both high-ionisation and 
high-density lead the high efficiency of gas-phase H2 formation 
and large D CI . 

By the way, we focus on the process that the ISM of dlrrs 
once is heated up by stellar-radiation, supernovae and/or other 
processes, and then cools down to the temperature at which star 
formatio n is effective u nder ISRF. Thus, re-calculating D a (eq. 
(35) of ldloverl d2003l) ) by using the physical parameters in 
which we interested, we have D a - > 10 _2 D M w at T > 5000 K 
as shown in the dotted line of figure ^ This critical value is 
on the same order of magnitude with t he typical valu e for dlrrs 
Ddtes ~ 1.3 X 10- 2 D MW dLisenfeld&Ferrarall998l). 

At T < 5000 K, the ISM cools off at a stroke. This is be- 
cause D CI becomes smaller than Aiim, and thus, H2 formation 
on dust-grain surfaces may dominate; therefore, yn 2 rapidly in- 
creases. Furthermore, if t he metallicity of d lrrs is higher than a 

-9 J ' II "I 

critical value ~ 10 ~Zq (Nishi & Tashiro 2000) (the same or- 
der with the metallicity of the most metal-deficient galaxy, I Zw 
18), metal cooling dominates H2 cooling at this temperature 
range. However, we are interested in the temperature at which 
the cooling timescale is the longest because global star forma- 
tion timescale is decided by the slowest physical process. Thus, 
the interesting temperature is found to be ~ 5000 K, at which 
gas-phase H2 formation dominates. Therefore, it is meaningful 
to study the SFHs of dlrrs on the assumption that gas-phase 
H2 formation is important. 

3. MODEL AND CALCULATION 

We calculate the time-dependent, non-chemical equilibrium, 
ionisation and molecule formation history of a gas system 
in dlrrs composed of hydrogen and helium. Neglecting the 
dynamics of the gas system in ISM, we concentrate on the 
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time evolution of the density and temperature in some re- 
gions in dlrrs. We consider the following 9 chemical species; 
H, H + , H , H 2 , H+, He, He + , He ++ , e . In dlrrs, the helium 
abundance is still maintained almost primordial value because 
dlrrs are still in an early stage of chemical evolution; therefore, 
it is assumed to be primordia l value (24% of hydrogen by mass 
fracti on) in our calculation ( Skil lman et alJll993t llzotov et all 
1999). 

At first, the ISRF in our model, J v , is assumed to have 
sam e energy dependence with that o f the solar vicinity, J® 
(see lMafhis. Mezger & Panagial ll983). and its intensity is nor- 
malised by that of the solar vicinity; J v — s • 7®. In order 
to investigate the dependence of yn 2 on the ISRF, we change 
the ratio s from 0.01 to 100 with the step of one order of 
magnitude. At next, decomposing the ISRF of solar vicinity 
into four parts, that is, early type stars, two types of disc stars 
(T eff = 7500 K, 4000 K) and red gian t stars (r eff = 3000 K) 
as in lMathis. Mezger & Panagial dl983l) . we investigate which 
degree each component affects vh, ■ Initial conditions of all cal- 
culations are Tm = l.Ox 10 6 K, « ini = 3.0 x 10 3 cm 3 , and the 
initial abundances of each chemical species are that in chemical 
equilibrium at 7^. 

4. Results 

Figure [2] shows vh 2 and ye- under some ISRF intensities. As 
the ISRF becomes more intense, these formative abundances 
monotonically decrease at a fixed temperature. vh 2 at T = 
5000 K, which is tightly connected with the SFRs of dlrrs in 
our model (see equation [6ji, reaches ~ 10~ 5 for the ISRF ten 
times more intense than that of the solar vicinity. This low 
H2 abundance means that the ISM cannot cool down and that 
star formation is delayed, and thus, this can explain the ob- 
served small SFRs of dlrrs in their quiescent phase. 

The dependences of yn, at 5000 K on the ISRF intensity 
are drawn in figure|3] The radiation from both disc stars (T e ff = 
7500 K, 4000 K) and giant stars (r eff = 3000 K) is as impor- 
tant as that from early type stars as the LW band photon sources 
in order to decrease vh 2 . The ISRF intensity, which is slightly 
stronger than that of solar vicinity, may be easily attainable 
because the ISRF intensity of the solar vic inity is still m uch 
weaker than that of the Magellanic Clouds ( Bl uhm & de BoeJ 
200lt lConfursi et alfeOOoV 

Here, we simply model the SFR, using Jeans mass, Mj, 
Jeans radius, rj, cooling time, f c00 i, and the mass conversion 
fraction from clouds to stars (star formation efficiency; SFE); 



SFR = 



SFE ■ Mj 

^cool 



In 2 . 



(5) 



This equation means that a cloud with mass M s and size r } will 
cool down in the timescale f c00 ] and some portions of the cloud 
which mass fraction are represented by SFE will convert into 
stars. Then, we apply the model to our simulation, and get 



SFR ~ l.Ox 10 



-4 (SFE\ 
10- 2 



[M /yr/kpc 2 ], (6) 




where vh 2 is the abundance of H2 at 5000 K, at which the cool- 
ing timescale is the longest. We can predict that the SFE of 



T[K\ 

Fig. 1. D CI as a function of the gas temperature T under some 
ISRF intensities are presented. The thick, long-dashed, short- 
dashed and dotted lines represent D CI for the ISRF intensity ra- 
tios s of 0.010, 0.10, 1.0, and 10.0, respectively. The ionisation 
degree and the density are self-consistently calculated. 

dlrrs is low (the order of 10~ 3 ), and ISRF of dlrrs is stronger 
than that of the solar vicinity in order to understand why the 
SFRs of dlrrs are such low. 



5. SUMMARY AND DISCUSSION 

We can conclude that the low mass stars, which are consid- 
ered to be H destroying photon sources, are as important as 
high mass OB stars as the sources of LW band photon in un- 
derstanding the episodic SFHs of dlrrs. That is, the radiation 
from many low mass stars burn in previous burst epoch of star 
formation destroys H , which is the intermediary of gas-phase 
H2 formation process, and suppresses further H2 formation 
with H2 dissociation, and keep the SFHs in quiescent phase. 
The schematic explanation of our conclusion is represented 
by figure Incidentally, we can say that low mass stars are 
more important than OB stars because the typical life time of 
OB stars is considerably short (~ 10 6 yr). Furthermore, low 
mass stars may be much more abunda nt in low s urface bright- 
ness galaxies than ordinary expected jLee et aljE( )04); that is, 
"bottom-heavy" initial mass function (IMF) is suggested there. 
If this IMF stands up in dlrrs, the contribution of the low mass 
stars to decrease y H , may become more and more important. 
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Fig. 2. yn 2 and vh- as a function of T under some ISRF inten- 
sities are presented. The four lines increased monotonically as 
the temperature decreases are vh 2 (labelled as "H2") and other 
four lines are y H - (also labelled as "H"). The thick, long- 
dashed, short-dashed and dotted lines represent the formative 
abundances for the ISRF intensity ratios s of 0.010, 0.10, 1.0, 
and 10.0, respectively. 
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the all components to be increased by an order of magnitude. 
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